The PhotoDissociation Region Toolbox: Software and Models for Astrophysical Analysis

The PhotoDissociation Region Toolbox provides comprehensive, easy-to-use, public software tools and models that enable an understanding of the interaction of the light of young, luminous, massive stars with the gas and dust in the Milky Way and in other galaxies. It consists of an open-source Python toolkit and photodissociation region models for analysis of infrared and millimeter/submillimeter line and continuum observations obtained by ground-based and sub-orbital telescopes, and astrophysics space missions. Photodissociation regions (PDRs) include all of the neutral gas in the ISM where far-ultraviolet photons dominate the chemistry and/or heating. In regions of massive star formation, PDRs are created at the boundaries between the H II regions and neutral molecular cloud, as photons with energies 6 eV $< h \nu <$ 13.6 eV photodissociate molecules and photoionize metals. The gas is heated by photo-electrons from small grains and large molecules and cools mostly through far-infrared fine-structure lines like [O I] and [C II]. The models are created from state-of-the art PDR codes that includes molecular freeze-out; recent collision, chemical, and photo rates; new chemical pathways, such as for oxygen chemistry; and allow for both clumpy and uniform media. The models predict the emergent intensities of many spectral lines and FIR continuum. The tools find the best-fit models to the observations and provide insights into the physical conditions and chemical makeup of the gas and dust. The PDR Toolbox enables novel analysis of data from telescopes such as ISO, Spitzer, Herschel, STO, SOFIA, SWAS, APEX, ALMA, and JWST.Comment: 22 pages, 10 figures, includes code listing

Dense, Parsec-Scale Clumps Near the Great Annihilator

We report on Combined Array for Research in Millimeter-Wave Astronomy and James Clerk Maxwell Telescope observations toward the Einstein source 1E 1740.7–2942, a low-mass X-ray binary commonly known as the "Great Annihilator." The Great Annihilator is known to be near a small, bright molecular cloud in a region largely devoid of emission in ^(12)CO surveys of the Galactic center. This region is of interest because it is interior to the dust lanes which may be the shock zones where atomic gas from the HI nuclear disk is converted into molecular gas. We find that the region is populated with a large number of dense (n ~ 10^5 cm^(–3)) regions of excited gas with small filling factors. The gas appears to have turbulent support and may be the result of sprays of material from collisions in the shock zone. We estimate that ~(1-3) × 10^5 M⊙ of shocked gas resides in our r ~ 3', Δv_(LSR) = 100 km s^(–1) field. If this gas has recently shocked and is interior to the inner Lindblad resonance of the dominant bar, it is in transit to the x_2 disk, suggesting that a significant amount of mass may be transported to the disk by a low filling factor population of molecular clouds with low surface brightness in larger surveys

Formation of Pillars at the Boundaries between H II Regions and Molecular Clouds

We investigate numerically the hydrodynamic instability of an ionization front (IF) accelerating into a molecular cloud, with imposed initial perturbations of different amplitudes. When the initial amplitude is small, the imposed perturbation is completely stabilized and does not grow. When the initial perturbation amplitude is large enough, roughly the ratio of the initial amplitude to wavelength is greater than 0.02, portions of the IF temporarily separate from the molecular cloud surface, locally decreasing the ablation pressure. This causes the appearance of a large, warm HI region and triggers nonlinear dynamics of the IF. The local difference of the ablation pressure and acceleration enhances the appearance and growth of a multimode perturbation. The stabilization usually seen at the IF in the linear regimes does not work due to the mismatch of the modes of the perturbations at the cloud surface and in density in HII region above the cloud surface. Molecular pillars are observed in the late stages of the large amplitude perturbation case. The velocity gradient in the pillars is in reasonably good agreement with that observed in the Eagle Nebula. The initial perturbation is imposed in three different ways: in density, in incident photon number flux, and in the surface shape. All cases show both stabilization for a small initial perturbation and large growth of the second harmonic by increasing amplitude of the initial perturbation above a critical value.Comment: 21 pages, 8 figures, accepted for publication in ApJ. high resolution figures available upon reques

Dense, Parsec-Scale Clumps near the Great Annihilator

We report on Combined Array for Research in Millimeter-Wave Astronomy (CARMA) and James Clerk Maxwell Telescope (JCMT) observations toward the Einstein source 1E 1740.7-2942, a LMXB commonly known as the "Great Annihilator." The Great Annihilator is known to be near a small, bright molecular cloud on the sky in a region largely devoid of emission in 12-CO surveys of the Galactic Center. The region is of interest because it is interior to the dust lanes which may be the shock zones where atomic gas from HI nuclear disk is converted into molecular gas. We find that the region is populated with a number of dense (n ~ 10^5 cm^-3) regions of excited gas with small filling factors, and estimate that up to 1-3 x 10^5 solar masses of gas can be seen in our maps. The detection suggests that a significant amount of mass is transported from the shock zones to the GC star-forming regions in the form of small, dense bundles.Comment: 26 pages, 7 figures, accepted for publication by the Astrophysical Journal, abstract abridge

A Resolved Ring of Debris Dust around the Solar Analog HD 107146

We present resolved images of the dust continuum emission from the debris disk around the young (80-200 Myr) solar-type star HD 107146 with CARMA at λ = 1.3 mm and the CSO at λ = 350 μ. Both images show that the dust emission extends over an approximately 10" diameter region. The high-resolution (3") CARMA image further reveals that the dust is distributed in a partial ring with significant decrease in a flux inward of 97 AU. Two prominent emission peaks appear within the ring separated by ~140° in the position angle. The morphology of the dust emission is suggestive of dust captured into a mean motion resonance, which would imply the presence of a planet at an orbital radius of ~45-75 AU

TADPOL: A 1.3 mm Survey of Dust Polarization in Star-forming Cores and Regions

We present {\lambda}1.3 mm CARMA observations of dust polarization toward 30 star-forming cores and 8 star-forming regions from the TADPOL survey. We show maps of all sources, and compare the ~2.5" resolution TADPOL maps with ~20" resolution polarization maps from single-dish submillimeter telescopes. Here we do not attempt to interpret the detailed B-field morphology of each object. Rather, we use average B-field orientations to derive conclusions in a statistical sense from the ensemble of sources, bearing in mind that these average orientations can be quite uncertain. We discuss three main findings: (1) A subset of the sources have consistent magnetic field (B-field) orientations between large (~20") and small (~2.5") scales. Those same sources also tend to have higher fractional polarizations than the sources with inconsistent large-to-small-scale fields. We interpret this to mean that in at least some cases B-fields play a role in regulating the infall of material all the way down to the ~1000 AU scales of protostellar envelopes. (2) Outflows appear to be randomly aligned with B-fields; although, in sources with low polarization fractions there is a hint that outflows are preferentially perpendicular to small-scale B-fields, which suggests that in these sources the fields have been wrapped up by envelope rotation. (3) Finally, even at ~2.5" resolution we see the so-called "polarization hole" effect, where the fractional polarization drops significantly near the total intensity peak. All data are publicly available in the electronic edition of this article.Comment: 53 pages, 37 figures -- main body (13 pp., 3 figures), source maps (32 pp., 34 figures), source descriptions (8 pp.). Accepted by the Astrophysical Journal Supplemen

Dynamically Driven Evolution of the Interstellar Medium in M51

Massive star formation occurs in giant molecular clouds (GMCs); an understanding of the evolution of GMCs is a prerequisite to develop theories of star formation and galaxy evolution. We report the highest-fidelity observations of the grand-design spiral galaxy M51 in carbon monoxide (CO) emission, revealing the evolution of GMCs vis-a-vis the large-scale galactic structure and dynamics. The most massive GMCs (giant molecular associations (GMAs)) are first assembled and then broken up as the gas flow through the spiral arms. The GMAs and their H_2 molecules are not fully dissociated into atomic gas as predicted in stellar feedback scenarios, but are fragmented into smaller GMCs upon leaving the spiral arms. The remnants of GMAs are detected as the chains of GMCs that emerge from the spiral arms into interarm regions. The kinematic shear within the spiral arms is sufficient to unbind the GMAs against self-gravity. We conclude that the evolution of GMCs is driven by large-scale galactic dynamics—their coagulation into GMAs is due to spiral arm streaming motions upon entering the arms, followed by fragmentation due to shear as they leave the arms on the downstream side. In M51, the majority of the gas remains molecular from arm entry through the interarm region and into the next spiral arm passage